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O ' Abstract 

o. 

^D . The melting parameter Tm of a classical one-component plasma is estimated 

using a relation between melting temperature, density, shear modulus, and crystal 
coordination number that follows from our model of dislocation-mediated melting. 
We obtain T^ = 172 it 35, in good agreement with the results of numerous Monte- 

'^ ! Carlo calculations. 
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1 Introduction 

The classical one-component plasma (OCP) is an idealized system of mobile ions of charge 
Ze, number density n, and temperature T, immersed in a neutralizing background of uni- 
form charge density —Zne. The OCP is realized in nature only at the enormous densities 
occuring in white dwarfs and neutron stars. The thermodynamics of the classical OCP is 
completely described in terms of the dimensionless coupling parameter ]l|] 
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where a = (3/47rn)^/^ is the Wigner-Seitz radius. In the quantum regime, one more 
parameter, a or T, is needed to characterize the system. Melting of a classical OCP occurs 
at a fixed value, F^, of the plasma coupling parameter. When F > F^, a OCP is either 
a glass 0, or it has a bcc crystal structure provided that it is subject to only hydrostatic 
stress. The evaluation of F^ for melting from the bcc structure has been the subject 
of extensive Monte Carlo (MC) calculations [3-16] employing the Ewald potential, which 
yields data pertinent to an infinite system from simulations using only a finite number 
of particles confined to a cubic computational cell with periodic boundary conditions. 
By fitting simple functional forms, guided by theory, to the measured excess potential 
energy per particle for both liquid and solid phases of the OCP, it is possible to obtain 
the Helmholtz free energy as a function of F. The intersection of the liquid and solid 
free-energy curves gives the value of the melting parameter F^. In their pioneering study 
PI, Brush, Sahlin and Teller observed melting in a 32-particle system at F^^ ^ 125. 
Subsequently Hansen Q] and Pollock and Hansen [^ followed with an improved calculation 
and found F^ = 155 ± 10. Van Horn |^ used the empirical Lindemann melting criterion 
to obtain F^ = 170 ± 10. Other MC studies resulted in the following values of F^ : 
144 0], 168 ± 4 §, 178 ± 1 §, 180 ± 1 0, 178 0, 172 |T2[, and 173 0. Values 



of very similar magnitude have been obtained in MC simulations of a strongly-coupled 
screened-Coulomb (Yukawa) system in the limit of zero screening: 171 |T^ and 171.8 |I5| . 



Recent path-integral MC simulations of the OCP |T6[ give Tm = 175. Hence, numerous 
MC studies suggest that F^ ~ 170 — 180 for the classical bcc OCP. 

In this paper we calculate F^ using a melting relation obtained from our model of 



dislocation-mediated melting [|1^, |18[. Before proceeding with the calculation of F^, we 
briefiy recapitulate the main ideas and assumptions of our melting model. As first pro- 
posed by Mott [|1^], dislocations are assumed to be the basic degrees of freedom underlying 



the melting transition. Dislocation interactions beyond a distance of order the mean dis- 
location separation are assumed negligible because of screening, and steric interactions 
are ignored. Accordingly, dislocations are taken to be non-interacting and therefore un- 
correlated, and are modeled as lines lying along the nearest-neighbor links of the lattice. 
The links coincide with the shortest perfect-dislocation Burgers vectors, which have mag- 
nitude b. The dislocation configurations (Brownian, self- avoiding, open, closed, etc.) are 
parametrized by a single parameter q > 1, in terms of which the mean dislocation length 
is given by (L) = Aqh/{q — 1). In addition to g, the partition function depends on the 
temperature-dependent effective dislocation line tension, that is, the energy cost to create 
unit length of dislocation at temperature T. The effective line tension vanishes at the 
critical temperature kBT^r = (Tb/ln{z — 1). Here z is the coordination number of the lat- 
tice and a, which we discuss in more detail below, is the p-dependent self-energy per unit 
length, p being the dislocation density. Dislocations proliferate as Tcr is approached from 
below, while at temperatures just above Tcr the partition function diverges, an indication 
that a new phase appears. So T^r corresponds to a phase transition, namely melting, 
and we identify T^r with the melting temperature, T^. A full defect theory of melting, 
a version of which is presently available |2^, would have to include the effects of both 
dislocations and disclinations. In our model we ignore the effects of disclinations under 
the assumption that they will produce only small changes of the order of 10% to the 



melting temperature. 

Under the assumption that dislocation strain fields are screened away at distances 
beyond the mean interdislocation spacing, the self-energy per unit length is given by |2T 
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where 2R ^ ^/ y/p is the mean distance between dislocations, z/ is the Poisson ratio, G 
is the shear modulus and a accounts for non-linear effects in the dislocation core. Many 



authors [^ have successfully used this ln(l/p) form for a, so we chose it as well, even 



though it has not been thoroughly investigated theoretically. Careful derivations |]23 
have been carried out only for nearly parallel dislocations. However, the ln(l/p) form is 
expected to hold in a three-dimensional ensemble of non-directed dislocations provided 
the mean dislocation length is much larger than the mean distance between dislocations, 
that is, (L)y/p ^ 1. In our model the ln(l/p) self-energy leads to a dislocation free energy 
F = —aip\np — a2p — asp""^, and the plnp term results in a first-order melting transition. 
We obtain the following melting relation: 
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(3) 



Here vws is the Wigner-Seitz volume, A = b^ /vws is a geometric factor characterizing 
the lattice, and piT^n) is the dislocation density at melt. Note that the factor ln(z — 1) 
explicitly accounts for the infiuence of crystal structure on melting. This melting relation 
plus experimental data on over half the elements in the periodic table gives b'^p{Tm) = 

0.61 ±0.20 ror 



In ref. |17| we applied Eq. (3) to the zero-pressure elemental data for more than half of 
the periodic table and found that it is accurate to 17%. Here we investigate the validity 
of this relation in the OCP limit by using it to calculate the value of F^, which is then 
compared to the available MC data. 

Calculation of F^ from Eq. (3) requires that we make the reasonable assumption that 
a^/6^p(Tm) is a pressure-independent constant. Then we can estimate this quantity for 
the OCP from zero-pressure data on the alkali metals. It is well known that the deviations 
of alkali-metal Fermi surfaces from perfect spheres are of order 1% or less, clear evidence 
that the valence electrons are very nearly free. In addition, the ratio of ionic radius to half 



the interatomic distance increases from 0.4 in Li to only 0.7 in Cs p4[, hence the overlap 
between alkali ions is small, and so to a good approximation the ions are effectively point 
charges. With respect to many of its physical properties (third-order elastic constants are 
one exception |^), an alkali metal can be regarded as a bcc lattice of point positive ions 
in a uniform background of free electrons, i.e., a one-component plasma. 

2 Analysis of alkali metal data 

Let us first discuss the temperature dependencies of G and vws, since their values to be 
used in Eq. (3) should be those at T = T^, not the measured values at room temperature. 



The fixed-pressure ratio of Wigner-Seitz volumes at T^ and T = is equal to 1 + PTm, 
where /3 is the volume expansivity. At j9 = 0, /5 is typically of order 10~^ K~^, and melting 
temperatures are at most about 4000 K, so vws changes by only a few percent between 
T = and T^. We can therefore always use room-temperature values for vws- 

In contrast to vws, the dependence of G on T is not necessarily weak. Its T-dependence 
involves two characteristic temperatures, namely the Debye temperature, T^, and the 
melting temperature. G is always monotonically decreasing with T, and is nonlinear for 
T ~ To and linear from To to T^. An accurate representation of G{T) at fixed density 
is achieved by ignoring the low-temperature non-linearity and approximating G{T) as a 
linear function of the reduced temperature T/Tm with the correct value G{0) at T = 

G{T) = G(0) (l - ^^) . (4) 



This straight-line representation turns out to be quite accurate: the maximum deviation 
of the data from the corresponding fitted lines is ~ 11% for the 22 metals analyzed in 

Hence, as follows from Eq. (3), 

1 _ 1 - Z//2 1 - 7 ^(O)^;^^^ .g. 
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where we have taken z = 8. For our analysis, we use Li, Na, K, Rb, and Cs, and omit 
Fr for which lattice constant data are not available. For the remaining 5 alkali metals, 
Vws = ct^/2, where the values of the lattice constant, a, are taken from [^. The values 



of both G(0) and T^ come from p8[. We take the values of 7 from ref. ||26[ for Na, K 



and Rb, and for Li and Cs we use 7 = 0.23, which is the average value over the 22 metals 



analyzed in [^. The corresponding values of u are taken from ref. |^. Averaging over 



the five alkali metals gives 

1 



0.385 ±0.052, (6) 



Aln(«V46V(Tm)) 
where the error is the root-mean-square deviation. 

3 OCP melting parameter F^ 

We first consider the Poisson ratio in the OCP limit. In terms of G and the bulk modulus, 
B, the Poisson ratio is given by |2^ 

W-2G 
"" 2{3B + G)' ^' 

We approximate B^'^^ by the bulk modulus of the electron gas since the negative electro- 
static (Madelung) contribution never exceeds 10% of the bulk modulus of the gas. The 
variation of 50CP -^j^j-^ density changes from a n^/^ dependence in the non-relativistic case 



to a n^/^ dependence in the extreme relativistic limit, whereas the OCP shear modulus 
always varies as n^/^. Hence, in a non-relativistic gas, B /G ~ n^/^ ^ 1, so i/ — >■ 1/2, as 
seen from Eq. (7). Although B and G both vary as n^^^ in the extreme relativistic limit, 
we find B/G ^ 1000 Z'^/^, so again B/G > 1 and i/ ^ 1/2. 

Hence, as follows from Eqs. (1),(3) and (6) with u = 1/2, the OCP melting parameter 
is given by 



(0.385 ± 0.052) Stt In 7 
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GOCP(r„) ' 

where we have used vws = l/''^- 

The bcc OCP elastic constants were recently obtained by Ogata and Ichimaru, us- 



ing MC simulations [^], as functions of P. However, the formula for the effective shear 
modulus used in ref. pQJ , 

Cll — Ci2 + 3C44 



G, 
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(9) 



is in fact the Voigt (upper) bound [^ 



on the shear modulus, and therefore does not give 
the correct value of G, which is known to always lie between the Voigt and the Reuss 
(lower) [^ bounds. 



An analysis by Kroner |3^ shows that successively narrower bounds can be placed on 
the shear modulus as the degree of disorder in grain orientation increases. In the limit of 
perfect disorder, the shear modulus can be obtained as a root of a cubic equation with 
coefficients that depend on the single- crystal elastic constants. In the case of the OCP, 
where the shear modulus is down by a factor of n^^^ from the bulk modulus, the cubic 
equation reduces to a quadratic with only one positive real root: 



G 



Cm + \Jcl^ + 12 (cii - C12) C44 



(10) 



In Table I we present the values of the elastic constants from ref. |^9| and the correct 
values of G°*-'^ as calculated from Eq. (10). 



p 


(cn - Ci2)/2 


C44 


qOCV- 


00 


0.02454 


0.1827 


0.0930 


800 


0.024(2) 


0.174(1) 


0.089(12) 


400 


0.025(2) 


0.167(1) 


0.087(11) 


300 


0.025(3) 


0.157(4) 


0.084(19) 


200 


0.019(3) 


0.12(1) 


0.064(28) 



Table I. The elastic constants and shear modulus, 
in units of (47r/3)i/3(Ze)2n4/3. 

Let us again assume a linear temperature dependence of G*-*^^ on T oc 1/P 

GOCP(r) = GOCP(oo)('l-^V 



(11) 



Fitting the values in Table I to this linear formula, and taking into account their uncer- 



tainties, we obtain |^| 



77 = 36.7 ±30.4. (12) 

Finally, we evaluate the OCP melting parameter F^ in the framework of melting as 
a dislocation-mediated phase transition. As follows from Eqs. (8), and (11), (12) with the 
value of G°'"^(oo) from Table I, 

F^ = 172±35. (13) 

This value is in good agreement with the available data from MC simulations, albeit 
with 20% uncertainty. We note that most (^ 2/3) of this uncertainty comes from the 
uncertainty in the value of r]. 

The OCP value of the parameter 7 defined in Eq. (4) is simply related to rj and F^ : 

V ^ 7°^"r^. (14) 

From Eqs. (12) and (13) we get 7^"-^^ = 0.21 ± 0.18, which is consistent with the value of 
7 at p ~ 0, namely 0.23 pef. 



4 Concluding remarks 

Our central value for F^, that is 172, agrees well with the more recent MC results. Two- 
thirds of the 20% uncertainty in this value is attributable to the error in the MC-calculated 
temperature dependence of the OCP single-crystal elastic constants. 

Our previous study of the melting curves of 18 elements [03 revealed that the melting 



relation (3) is in good agreement with data up to pressures ~ 100 — 200 GPa. Here we have 
demonstrated that Eq. (3) also holds in a classical OCP. These successful comparisons 
of Eq. (3) with experimental data and MC calculations suggest, but of course do not by 
themselves prove, that melting is a dislocation-mediated phase transition. 
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